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Abstract 

We present a detailed study of the interplay between chiral sym- 
metry and spectral properties of the Dirac operator in lattice gauge 
theories. We consider, in the framework of the Schwinger model, the 
fixed point action and a fermion action recently proposed by Neu- 
berger. Both actions show the remnant of chiral symmetry on the 
lattice as formulated in the Ginsparg- Wilson relation. We check this 
issue for practical implementations, also evaluating the fermion con- 
densate in a finite volume by a subtraction procedure. Moreover, we 
investigate the distribution of the eigenvalues of a properly defined 
anti-hermitian lattice Dirac operator, studying the statistical proper- 
ties at the low lying edge of the spectrum. The comparison with the 
predictions of chiral Random Matrix Theory enables us to obtain an 
estimate of the infinite volume fermion condensate. 
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1 Introduction 



A local lattice theory describing the physical degrees of freedom of the corre- 
sponding continuum theory necessarily breaks chiral symmetry £Q. In a usual 
discretization, like the Wilson action, this breaking is so violent, that no trace 
of the chiral properties of the continuum is kept in the lattice theory This 
involves many inconveniences - in particular for lattice QCD. At the classical 
level, for fixed gauge field configurations, zero modes of the (Euclidean) Dirac 
operator have no definite chirality, and the Atiyah-Singer theorem of the con- 
tinuum theory j2] , relating the index of the Dirac operator to the topological 
charge of the background gauge configuration, finds no strict correspondence 
on the lattice. At the quantum level, the chirality-breaking terms - when 
combined with the ultraviolet divergences of the theory - give raise to non- 
universal finite renormalizations. Examples are the additive renormalization 
of the bare quark mass m q and the finite multiplicative renormalizations of 
the chiral currents, which spoil the usual current algebra. Since the chirality 
is explicitly broken, the definition of an order parameter for the spontaneous 
breaking of the chiral symmetry is not straightforward. 

In an early paper [3| of lattice quantum field theory, Ginsparg and Wilson 
made definite the concept of a chiral limit in the framework of a lattice theory 
breaking explicitly chiral symmetry, providing a general condition - we will 
refer to it as the Ginsparg- Wilson Condition (GWC) - for the fermion matrix 
of the lattice theory, i.e. the lattice Dirac operator. For a long time Ginsparg 
and Wilson's ideas remained academic, since no acceptable solution of the 
GWC was found in the case of dynamical gauge fields, locality of the lattice 
action being the real bottle-neck in this matter. 

New attention on the GWC was raised by Hasenfratz 4], who pointed out 
that the fixed point (FP) action, which is local by construction, satisfies this 
relation. This result is natural, since the FP action is a classically perfect 
action. 

Chiral symmetry on the lattice is made explicit in the overlap formalism 
recently Neuberger [B] found a form of the corresponding Dirac operator, 
which also satisfies the GWC |Zj. 

In a series of papers, the GWC was theoretically analyzed, showing that 
it is a sufficient condition for the restoration of the main features of the 
continuum (symmetric) theory. In the following we refer to actions/fermions 
satisfying the GWC as to GW actions/fermions. At the classical level, for FP 
actions, the Atiyah-Singer theorem finds correspondence on the lattice [HE]; 
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at the quantum level, no fine tuning, mixing and current renormalization 
occur, and a natural definition for an order parameter of the spontaneous 
breaking of the chiral symmetry is possible [5J. The explicit form of the 
lattice symmetry corresponding to the chiral symmetry of the continuum 
has been identified On the basis of this symmetry it has been shown 

[TT] that the low-energy mechanisms of the continuum QCD (e.g. Goldstone's 
theorem for pions, solution of the U(l) problem) emerge also in the theory 
with finite cut-off. 

Monte Carlo calculations require actions with finite number of couplings 
(ultra-local actions), while any solution of the GWC - even when local - 
is expected JT2J to be extended all over the lattice. In the case of the FP 
action, the renormalization group theory ensures exponential damping of the 
couplings with the distance, and a parametrization procedure is viable in 
principle; the point to be verified is to what extent an approximation of the 
FP action in terms of a restricted set of couplings is able to reproduce the 
nice theoretical properties of the 'ideal' FP action. 

In this paper we check these issues in the case of the Schwinger model 
(d = 2 QED for massless fermions), testing a parametrization of the FP 
action for the non-overlapping block-spin transformation. We focus mainly 
on the interplay between chiral symmetry and spectral properties of the Dirac 
operator. Preliminary results were presented in |2] and ^Hj. We also test 
Neuberger's proposal [0]. Previous investigations on this Dirac operator in 
the framework of the Schwinger model were accomplished in ^Hj and (for 
an application of the standard overlap formalism, see JH])- Neuberger's Dirac 
operator is determined through a configuration-wise operatorial (numerical) 
projection [7j of the Wilson operator. In this case an explicit parametrization 
of the Dirac operator is not yet available; also the precise locality properties 
of this operator have not been established (see, however |19j). 

Based on universality argument, suggested by Leutwyler and Smilga's [20] 
sum rules for the eigenvalues of the chiral Dirac operator in a finite volume, 
Shuryak and Verbaarschot j.21J hypothesized that statistical properties of the 
spectrum of this operator at its lower edge are described by Random Matrix 
Theory (RMT, for a recent review cf. [22] )■ This ansatz has been explicitly 
verified [23J. When chiral symmetry holds in the lattice theory, the chiral 
version of RMT (chRMT) should apply [2H [21 123- So, another way to 
check the effective restoration of the chiral symmetry by GW actions, is to 
compare [2B] the statistical properties of their spectrum with the prediction 
of chRMT. In this paper we also address this point. 
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The sequel of the paper is organized as follows. First we summarize the 
main properties of GW lattice fermions in relation to chiral symmetry; in 
Section |3] we report our numerical results about the spectrum of the Dirac 
operator, in particular the chirality of (quasi) zero modes, their dispersion 
on the real axis, and agreement of the index of the Dirac operator with the 
geometric definition for the topological charge. We also show the results of 
the calculation of the chiral condensate in a finite volume according to the 
prescription given in jS] . In Section 0] we concentrate on the study of the 
statistical properties of the eigenvalues of the Dirac operator, investigating 
the universality class at the lower edge of the spectrum. Comparison with 
the predictions of chRMT allows us to extract the infinite volume fermion 
condensate from lattice data. The last section is devoted to discussion and 
conclusions. 



2 Ginsparg- Wilson fermions 

The GWC for the fermion matrix T> xx i describing massless fermions reads 

\{V X , X ,^} = (Vj 5 RV) XtXl , (1) 

where R x , x > is a local matrix in coordinate space, i.e. whose matrix elements 
vanish exponentially with the distance. Usual lattice actions obey the general 
7 5 -hermiticity property = 7 5 X>7 5 . In the particular case R x , x ' — \&x,x' 
the GWC then assumes the elegant form 

V + £> f = V ] V =VV ] . (2) 

Fixed point action. Any FP action of a block spin transformation (BST) 
satisfies the GWC [3] . Such an action is the FP of a recursion relation of the 
form, 

-kI L (m 1 r-) "lyiU) , (3) 

where U(V) is fixed by the solution of the pure-gauge part of the FP prob- 
lem, uj(U) is the matrix in space-time and color indices defining the (gauge- 
covariant) block-spin average, and free parameter. The operator R XX ' 
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has an expression [3] in terms of ui(U), and it can in general be assumed to 
be trivial in Dirac space. In the case of the non-overlapping BST considered 
here and in the particular form (J2J) is fulfilled. 

Neuberger's Dirac operator. In this approach jHj the starting point is 
the Wilson Dirac operator at some value of m G (—1,0) corresponding to 
2^ < k < 2D-2 ( K * s ^ ne well-known hopping parameter, not to be con- 
fused with the parameter of the block-spin transformation) and then one 
constructs 

V Ne = 1 + 75 e( 75 V m ) , where e(M) = -J£= . (4) 

y M 1 

From the definition (J3J) it is evident that V^ e satisfies (j2j) 0. Here we 
consider the case k = |, m = — 1 (for a discussion in the context of this 
study cf. HE]). 

2.1 Chiral properties 

Dirac operators T> obeying the GWC do have non-trivial spectral properties. 
For simplicity we consider here operators, like Ppp and T>^ e , satisfying the 
special version (J2J of the GWC. 

i. [V , Z>t] = 0, i.e. V is a normal operator; as a consequence, its eigen- 
vectors form a complete orthonormal set. 

ii. The spectrum lies on a unit circle in the complex eigenvalue-plane cen- 
tered at A = 1. 

iii. The property (i), together with the hermiticity property of the fermion 
matrix, implies: 

7 5 "a oc , if A ^ A , 

7 5 w A = ±v x , if A = AGR. (5) 

Here v\ denotes an eigenvector of T> with eigenvalue A. Just as in 
the continuum, the eigenvectors of complex-conjugate eigenvalues form 
doublets related through 7 5 ; moreover all real-modes have definite chi- 
rality. For the general form (pQ) of the GWC this property holds only 
for the zero modes. 
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2.2 The At iyah- Singer theorem on the lattice. 

The Atiyah-Singer index theorem (ASIT) in the continuum relates the topo- 
logical charge Q(A) of a differentiable gauge field configuration A to the 
difference between the numbers of positive and negative chirality zero modes 
of the Dirac operator, 

Q(A) = mdex(A) = n + - n_ . (6) 

On the lattice, an index for T> may be defined in a way analogous to the 
continuum, explicitly expressed by the relation [Hj 

index(tf) = -tr( 7 5 RV (U)) . (7) 

In the case of GW Dirac operators satisfying the special GW condition (j2J, 
the above relation comes out trivially considering that (due to their proper- 
ties listed in the previous section) only the modes with real non-vanishing 
eigenvalues contribute to the trace in the r.h.s, where R = 1/2; since the 
overall chirality must be zero, it reproduces up to a sign index(£7). 

This index([7) can be used to define a fermionic lattice topological charge 

Q ferm ([/) = index(f/) (8) 

for which the ASIT is satisfied by definition. 

In the case of the FP action the fermionic definition (|SJ) coincides jH] 
with the pure-gauge quantity Qp p (U), the FP topological charge [2B] of the 
configuration U: 

Q lcvm (U) = index(f/) = Q Fp (U) . (9) 

The non-obviousness of this relation relies on the fact that Qf p (U) can be 
defined a priori in the pure gauge theory, without any regard to the fermion 
part. We stress that this result is particular for a FP action, having no 
counterpart for a general (non-FP) GW action. Of course, in practical im- 
plementations one relies on approximate parametrizations of the FP Dirac 
operator and the strictness of the relation is consequently lost. 

Not all fermionic formulations of the topological charge, associated to 
different lattice Dirac operators, are equivalent, since they can be more or 
less reminiscent of the continuum 29J. The FP definition is in this sense the 
most reliable, since it is essentially based on a renormalization group guided 
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procedure of interpolation (but maybe difficult to implement in realistic sit- 
uations like QCD). Here, as a general criterion, we compare the results for 
T>p p and T>^ e with the geometric definition 

Q geo (U) = ± J2 fc*ln(M*)) > ( 10 ) 

X 

which is in this context the most natural one and the closest to the continuum 
definition. In d = 2 the geometric definition is the FP for a particular BST 

m 

In two dimensions, a theorem of the continuum - the so-called Vanishing 
Theorem [~~0] - ensures that only either positive or negative chirality zero 
modes occur. It seems plausible that this theorem applies for FP actions 
[2*7] : we know of no proof that it should apply to GW fermions in general. 



2.3 The fermion condensate 

In |3j Hasenfratz proposed a subtraction procedure for the fermion conden- 
sate inspired by (JTJ) (an analogous prescription for a physical lattice fermion 
condensate was proposed by Neuberger in an earlier paper |3T| in the specific 
framework of the overlap formalism and for the particular case R = 1/2). 
The definition of the subtracted lattice condensate reads 

(« S ub = -4 (to(^-R)) » (11) 

/ / gauge 



V \ \v 



where V is the (finite) space-time volume. The expectation value on the right- 
hand side denotes the gauge averaging including the determinant weight, 
corresponding to the full consideration of dynamical fermions. With (jTJ we 
may rewrite the right-hand side as 

gauge 



2V \ \V Z>t 



Because of 75-hermiticity the trace in the above expression vanishes, except 
when a zero mode of T> occurs, in which case a regulator-mass fi must be 
introduced, 

V -> T)(ji)=V +fil. (13) 

In d = 4 for more than one flavor, Nf > 1, (ip^sub is an order parameter 
for the spontaneous breaking of the chiral symmetry: the contribution of the 
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zero modes to the gluon average vanishes when \x — > because of the damp- 
ing effect of the fermion determinant. We conclude that for finite volume 
(V ; V ; )sub(A i ) vanishes in the chiral limit [3]. 

Here we consider d = 2 and Nf = 1; the chiral symmetry is broken by 
the U(l) anomaly Still we want to use the above expressions to study the 
condensate. In this situation (ipip) SVi b ^ even in a finite volume. The 
configurations responsible for the non-zero fermion condensate in a finite 
volume are those from the \Q\ = 1 sector; indeed, if T> has just one zero 
mode - which is possible due to the ASIT and Vanishing Theorem only for 
\Q\ = 1 - the quantity (detZ> trD -1 ) has a non-zero limit when /j, — > 0. 
The effect of the subtraction is in general (for any Nf) just to remove the 
contribution of the Q = sector. 

The situation in the infinite volume, when the fermion condensate is 
obtained through the sequence of limits lim M ^ nm v^oo ; is different. In this 
case the role of the exact zero modes is irrelevant, the quasi zero modes 
(the modes with small but non- vanishing eigenvalues) are responsible for the 
non-zero fermion condensate as the Banks-Casher formula [32] shows. 

We observe (cf. [HIE! as well as the discussion in |3S|) that (fTTj) may be 
rewritten in the form 

(■H>)*xb = ~ (t r l i J ) with V =V{l-RV)- 1 . (14) 

gauge 



v \ \v 



The redefined fermion matrix T> is anti-hermitian and anti-commutes with 
f: 

T> ] = -V , {P, 7 5 } = 0, (15) 
from which it follows that T> has a purely imaginary spectrum, 

Vv = -i\v, AGl. (16) 

The zero modes of T> coincide of course with those of T> . The replacement 
T> — > T> makes manifest the chiral invariance implicit in the original operator 
(observe however that according to Hasenfratz's subtraction prescription the 
dynamics of fermions is still given by the chirality-breaking matrix T>). 
The spectral density of T> , 

/ > \ 1 dN 
" (A) = V (17) 
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is obtained from the eigenvalue spectrum sampled with the weight of (detT> ) Nf 
It complies with the Banks-Casher formula for the (infinite-volume) sub- 
tracted fermion condensate 



lim lim (^^) SU b(A t ) = h m h m p(A) ; (18) 

the sequence of the limits is essential. 

In the case of the simpler GWC (J2J), the spectrum of V is obtained by 
mapping the spectrum of T> (which in this case, we recall, lies on a unit 
circle in the complex plane) onto the imaginary axis by the stereographic 
projection: 

y ' 



x - M 1 - 2 ) • ( 19 ) 



The spectral density agrees with that of T> in 0(A). 



2.4 Spectral microscopic fluctuations 

In the framework of QCD with arbitrary Nf, Leutwyler and Smilga [20] 
derived, through the calculation of low-energy effective lagrangians, sum 
rules which are obtained in a limit where QCD reduces to a simple matrix 
model. These arguments suggest [21] , that the symmetry of the theory (in the 
simplest case just the chiral one) rather than its detailed dynamical structure 
dictates the microscopic fluctuations of the spectrum of the massless Dirac 
operator; they should have a universal behavior in a scaling variable z = VEX 
in the limit V — *■ oo. Here E denotes a free parameter related to the dynamics 
of the system. For fixed z, X — > in the thermodynamic limit and the part of 
the spectrum near zero - which is also the most relevant for the continuum 
limit - is under study. 

The simplest model encompassing chiral symmetry is p!lj a chiral Random 
Matrix Theory. Since GW actions are expected to describe massless quarks 
and implicitly realize chiral symmetry on the lattice (explicitly for the matrix 
T> which anti-commutes with 75), it is natural to investigate whether the 
predictions of chRMT for the microscopic fluctuations of the spectrum apply; 
for GW fermions cf. [211 EH]- Previous studies for staggered lattice fermions, 
which partially restore chirality, where already accomplished for the SU(2) 
[S3 HUE EZl EH] and the SU(3) [US] lattice gauge theory in four dimensions. 

Three classes of universality are predicted within chRMT j2H], corre- 
sponding to orthogonal, unitary or symplectic ensembles (chOE, chUE or 
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chSE respectively); universality of related correlation functions has been 
proved in [30! • Starting from Leutwyler and Smilga's effective lagrangians, 
it is possible to show f22| that the microscopic fluctuations of QCD (with 
three colors) are described by chUE. In our context, QED 2 , no additional 
symmetry for the Dirac operator beyond the chiral one is present; in this 
situation one also expects chUE universality [2*5] . 

Spectral fluctuations may be studied through a variety of probes (for an 
incomplete list of such studies within lattice gauge theory, cf. jSH HUH EH EM 
). Here we consider just three of them: 

• the probability distribution of the smallest eigenvalue P(A min ); 

• the microscopic spectral density p s (z) 

p a {z)= lim Ip(^) , (20) 



where p(X) is the spectral density (fTTjl : 
• the number varianc^E 

E 2 (S , S) = ((N(S , S) - N(S , S)) 2 ) gauge , (21) 

i.e. the variance of the number of eigenvalues in the interval [So, S]. 

Different numbers of flavors Nf correspond to different predictions from 
chRMT. In the terminology of lattice calculations, Nf = corresponds to the 
quenched situation, while Nf = 1 (or more) to the dynamical one. There- 
fore comparison with chRMT allows us to check whether the dynamics of 
the fermions has been effectively included in the simulations I n the 

dynamical setup, the Banks-Casher relation ()18j) relates the parameter E to 
the (subtracted) fermion condensate in the infinite volume: 

E = — lim lim (^) sub . (22) 

II— >0 V^rOO 

In chRMT there are specific predictions concerning the spectral distribu- 
tion in different topological sectors. This theory knows nothing about the 
dynamical content of gauge fields; the parameter z/, which defines the topo- 
logical sector, is the difference between the numbers of left- and right-handed 



1 We thank J.-Z. Ma for drawing our attention on this particular test. 



10 



modes in the matrix representation. In the framework of gauge theories v 
corresponds to the index of the Dirac operator: for GW fermions it can be 
defined as in the continuum and produces a fermionic definition of the topo- 
logical charge on the lattice (see Section I2~T1 eq. (J7J) and (jHJ)). This is not 
the case for other (non-GW) lattice actions, where the index - and as a con- 
sequence the topological charge - is not well defined causing ambiguities in 
the comparison with chRMT [3%] . 

In our discussion of the numerical results we rely on dimensionless quan- 
tities. The lattice spacing a and the (physical) gauge coupling constant e are 
dimensionful and we assume the usual asymptotic relation to the dimension- 
less coupling P, 

= — . (23) 
e a 

In order to compare with the (theoretical) continuum values we therefore 
have to include corresponding factors of e. The physical lattice size (L a) is 
then in dimensionless units 

Lae = ^. (24) 
For Nf = 1 the Schwinger mass is 

1 

Ct 771 Schwinger / ■ v / 

VP 7T 

The continuum, infinite volume value for the condensate is 

-W)oant = ce where c = ^f P ^ sa 0.15989 , (26) 

2 7T V7T 

(7 denotes the Euler constant) and its dimensionless lattice partner (corre- 
sponding to (V^)cont/e) is 

W)lat ,7 /x /q M ^ 

ea = (W/lat VP ■ (27) 



3 The numerical analysis 

The lattice gauge fields U x ^ are in the compact representation with the usual 
Wilson plaquette action. 
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In the FP Dirac operator was parametrized as 

3 

v Fp (x, y) = J2J2 °i u ( x > f) > with y = x + 5f. (28) 

i=0 x,f 

Here / denotes a closed loop through x or a path from the lattice site x to 
y = x + Sf (distance vector 5f) and U(x, f) is the parallel transporter along 
this path. The aj-matrices denote the Pauli matrices for i = 1,2,3 and the 
unit matrix for i — 0. Geometrically the couplings have been restricted to lie 
within a 7 x 7 lattice. The action obeys the usual symmetries; altogether it 
has 429 terms per site. The action was determined as the (numerically ap- 
proximate) FP of the Dirac operator for gauge fields distributed according to 
the non-compact formulation with the Gaussian measure. Excellent scaling 
properties, rotational invariance and continuum-like dispersion relations were 
observed at various values of the gauge coupling (3 [T3|. Here we study the 
action (|28|) only for the compact gauge field distributions (see, however, the 
results in [Ej). In this case the action is not expected to exactly reproduce 
the FP of the corresponding BST, but nevertheless it is still a solution of the 
GWC; violations of the GWC are instead introduced by the parametrization 
procedure, which cuts off the less local couplings. 

The situation is different for Neuberger's operator, which is determined 
for each configuration on the basis of (J1J): in our simple context computer 
time was not really an obstacle and therefore we computed e(7s ZVi ) through 
diagonalization with machine accuracy (cf. ^Ej for more details, as well as 
[4*H 14*2*1 IT9"1 H3j for more efficient approaches in d — 4) . 

Essentially uncorrelated gauge configurations have been generated in the 
quenched setup. As a measure we have used the autocorrelation length for 
the geometric topological charge, which - for Metropolis updating - behaves 
like exp(1.67/3 — 3) for the range of gauge couplings studied. Unquenching 
is obtained by including the fermionic determinant in the observables; bet- 
ter ways to include the dynamics of fermions with the present actions - in 
particular Neuberger's operator - are still being developed |44j . From earlier 
experience [121 Ej we learned that our procedure is justifiable at least for the 
discussed model and the presented statistics. We perform our investigation 
on sets of 5000 configurations (for (3 = 2, 4) and 10000 configurations (for 
(3 = 6) for lattice size L = 16, and 5000 configurations (for (3 =4, 6) and size 
L = 24. 
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L -0.5 0.5 1.5 2.5 L -0.5 0.5 1.5 2.5 L -0.5 0.5 1.5 2.5 

Re 

Figure 1: Eigenvalues of the parametrized FP Dirac operator at the values 
of ft — 2,4, 6 (from left to right), on 16 2 -lattices and sampled over 25 gauge 
configurations each. 

3.1 The spectrum 

Fig. E] shows the spectrum of the studied FP Dirac operator collectively for a 
subset of 25 uncorrelated gauge configurations in thermal equilibrium. While 
Neuberger's Dirac operator has (within computational precision, in our case 
14 digits) exactly circular spectrum, in the case of the FP action the effect 
of the truncation in coupling space due to the ultra-local parametrization is 
displayed by an increased fuzziness for lower (3. 

The spectrum of the Dirac operator is gauge invariant and so it can be 
expanded in series of loop operators. At lowest order in the lattice spacing 
only few of them are independent, but when considering higher order correc- 
tions operators of larger (geometric) extent have to be introduced in order to 
have a complete set. So it may be argued that the truncation of the couplings 
to a 7 x 7 lattice for the present parametrization of the FP action implies an 
error in the form of an operator of some (high) dimension: 

App(ET) - A par (C7) = {k \U) , k = dim[0«] . (29) 

In order to estimate quantitatively the scaling behavior of the deviations 
of the spectrum from the ideal circular shape, we defined a mean deviation 
| A — 1| from the unit circle in an angular window of |arg(l — A)| < tt/A. The 
displayed behavior with f3 of the average width (standard deviation) a of 
that distribution is a oc 1/ ' (3 2A1 ~ a 5 , thus suggesting a dimension-5 operator 

Another effect observed for the FP action and due to the truncation, is 



13 




1.0 



0.5 



0.0 



1.0 



0.5 



0.0 



-0.50 -0.25 0.00 0.25 0.50 -0.10 -0.05 0.00 0.05 0.10 -0.02 -0.01 0.00 0.01 0.02 
zero mode position 



Figure 2: Distribution of the (real) eigenvalues for the quasi zero modes 
for the FP action on a 16 2 lattice; from left to right: f3 =2, 4, 6 (observe 
the different scales of the abscissa). For convenience the entries have been 
normalized to the maximum entry. 
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Figure 3: Distribution of the chirality (^,75^) of the quasi zero modes 
of the Dirac operator for the FP action on a 16 2 lattice; from left to right: 
/3 — 2, 4, 6 (observe the different scales of the abscissa). We show only 
the positive chirality contribution. For convenience the entries have been 
normalized to the maximum entry. 



the scattering of the zero modes on the real axis. The relative spread is tiny 
at large values of f3 but becomes relevant when decreasing [3. It damages the 
nice theoretical properties of the operator causing difficulties in the numerical 
approach (see the following). From Fig. El we see that the spread is ~ 0.002 
for /3 — 6, ~ 0.01 for j3 — 4, while for (3 = 2 it becomes substantial. 

A similar effect can be observed in the distribution of the chirality of 
the (almost) zero modes, displayed in Fig. El In this case, the accuracy of 



14 



Op. 


p(P = 2) 


piP = 4) 


p((3 = 6) 




74.22 


99.70 


100.00 


£>Fp 


96.58 


100.00 


100.00 




100.00 


100.00 


100.00 




91.00 


99.84 


99.96 



Table 1: Upper part: percentage p(/3) of configurations where the index of 
the Dirac operator index (U) agrees with the geometric charge. Lower part: 
percentage p(/3) of configurations where the Vanishing Theorem is fulfilled. 

realization of chiral symmetry is directly probed. For Neuberger's operator 
the chirality of the zero modes is ±1 within numerical precision. 

Table Ogives the results of the comparison between the index of the Dirac 
operator index(C7) (the modes are counted according to the sign of their 
chiralities) and the geometric charge of the gauge configuration. Observe 
that in the case of the FP action the agreement is excellent already at (3 = 2, 
thus indicating a close correspondence between the geometric definition and 
the FP topological charge. For this value of f3 the deviation is comparatively 
large for Neuberger's operator, reflecting in this respect features of the 
Wilson operator entering (jlj). 

The lower part of the table demonstrates that for T>^ e zeio modes have 
just one chirality here. For "Pppthis situation is also rapidly approached 
towards the continuum limit, as predicted by the Vanishing Theorem. 

3.2 The fermion condensate in a finite volume 

The finite-volume fermion condensate is determined following the prescrip- 
tions of Section 12.31 In order to check the correctness of our determination 
we verify the Ward identity (see also [TT]). valid in the limit \x — > 

-lim(#) S ub(/<) = limx'O) ( 30 ) 

With *V) = ^<( indeX ( f/ )) 2 > gaug e^)- ( 31 ) 

In Table El we compare the values for both sides of (J3Uj) (extrapolated 
at \i = 0): excellent agreement between the two quantities is displayed by 
data for Neuberger's action. In the case of the FP action, x'(/-0 i s unstable 
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Figure 4: The subtracted chiral condensate as a function of /i for Neuberger's 
(left) and FP (right) operator on a 16 2 lattice: (3 — 2 (circles), (3 = 4 
(diamonds), /3 = 6 (triangles); the full symbols denote the zero-mass estimate 
obtained through a linear extrapolation. For (3 > 4 they agree within the 
error bars with the finite volume numbers (indicated by the pointers on the 
left-hand side) from theory [4*5] . 
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Table 2: Subtracted condensate (^V)sub an d x! f° r the FP and Neuberger's 
action (all obtained in the limit \i — > 0). The last column indicates the 
continuum value in the corresponding physical volume |45j . 



for \x — > 0: this is simply explained by the spread of the zero modes which 
prevents the cancellation between the denominator factor /x with the lowest 
eigenvalue in the fermion determinant (which is \x for an exact zero mode). 
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The subtracted condensate for the (approximate) FP action is plagued by 
similar problems. For (3 = 4 and L = 24 the determination of the condensate 
has large error bars and for (3 = 2 and L = 16 it is not possible at all due 
to the instability of the \x — > limit. Except for these two cases the relative 
accuracy of our results is between 3 and 6 %. In Table El we also compare 
our data with the continuum value in the corresponding physical volume, 
assuming the asymptotic relation (j2HJ); the discrepancy is always smaller than 
one standard deviation. 

In Fig. H]we report {ipip) SU b(/-0 f° r a range of // values in the case of Neu- 
berger (left) and FP (right) action; the expected linear asymptotic behavior 
allows for the extrapolation at /i = 0. 

4 Statistical properties of the spectrum 

GW actions realize chiral symmetry on the lattice. This realization is im- 
plicit, in the sense that the Dirac operator breaks the chiral symmetry, as 
usually defined, by a local term 0(a). We make the symmetry manifest and 
construct a fermion matrix T> which anti-commutes with 7 5 . This is achieved 
through the transformation proposed in ()14j) . In this section we study the 
microscopic fluctuations of the spectrum of this operator. 

4.1 The trivial sector: Nf = and v = 

Since chRMT gives specific predictions for given topological number v and 
number of flavors Nf, it allows to disentangle the statistical properties of the 
spectrum of the lattice Dirac operators from the more subtle questions of the 
correct reproduction of the dynamics of fermions and of the identification of 
the topological charge with the index (see the discussion on this latter point 
in Section |23J). Consequently, we start our investigation in the simplest case 
N f = and u = 0. 

We expect that a gauge theory of the type studied here - without further 
symmetries - should lead to a spectral distribution in the universality class 
of chUE This is to be checked with the data. 

Comparison of the statistics of the Dirac spectrum distribution with 
chRMT requires the knowledge of the inherent scale parameter E that is 
related to the dynamical properties of the system. Technically S may be 
determined by comparing spectral data with chRMT predictions for a par- 
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Figure 5: The ratio between the xi (x 2 with chUE) and xl (x 2 with chSE) 
in the case of the FP action as a function of the physical lattice size L/y/]3. 

ticular statistics. The consistency of the determination can be then checked 
by using the obtained value of S as an input parameter for further compar- 
isons with chRMT. 

Smallest eigenvalue. A simple way to determine £ is to look at the proba- 
bility distribution of the smallest eigenvalue P(A m i n ). This statistics contains 
spectral information on a region of small enough values of A not to be affected 
by the macroscopic, non-universal (i.e. not described by chRMT) component 
of the spectrum. Of course the macroscopic component can be easily elimi- 
nated from lattice data through an unfolding procedure, which - at this level 
- would introduce unnecessary arbitrariness, however. 

As discussed, we expect results lying in the universality class of chUE. 
However, in order to decide on the applicability of chRMT in an unbiased 
approach, we tested the predictions jlH] for the distribution of the smallest 
eigenvalue from all three variants of chRMT through a standard best-fit 
procedure on our lattice data. This test clearly ruled out the chOE, which 
has a completely different distribution shape. The two remaining ensembles 
chUE and chSE seemed to fit data in disconnected - complementary - regions 
of the lattice parameter space (L,j3). 

In Fig. El we plot the ratio of x 2 -values obtained for the fits to the two 
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distributions as a function of the physical size. We find that in the region of 
large physical volumes chUE is preferable (the ratio is smaller than unity), 
while in the region of small physical volumes chSE provides better fits to 
the data. The 'transition' seems to take place at L/y/]3 ~ 9. In Fig. El we 
compare our lattice results for P(A m j n ) with the best fit prediction from chUE 
and chSE for two cases in the large-volume region (upper part) and in the 
small- volume one (lower part). In Table El we give the resulting values of S. 
Within chRMT this value is intrinsically volume-independent. The obtained 
value of £ is taken as an input parameter for the subsequent investigation. 

For small physical volumes we find chSE distribution shapes; however, 
this result is in conflict with other properties of such ensembles. For chSE 
the eigenvalues should come in degenerate pairs, which is not true in our 
cas^l 

Such a volume dependence of our results does not come unexpected and 
is related to the limitations of applicability of RMT. Actually, in QCD one 
expects RMT to be applicable [201 El] for 



From the continuum results we expect on the lattice ([25]). and then reads 
yfit <C Lfy/~j3, to be compared with the discussed observations. We suspect 
that the observed behavior at small volume is outside the scope of RMT and 
that one may be misled by the superficial agreement with chSE. We return 
to a discussion in the last section. 

Microscopic spectral density. Using £ as input parameter, as discussed, 
we checked the universal behavior of the microscopic spectral density p s {z). 
In this case the macroscopic unfolding of the lattice data was necessary (in 
this we followed essentially the suggestions in jHHl)- We report the results in 
Fig. [7J again comparing with the chUE [2I| and chSE [171 EH] predictions. 
We see a (universal) behavior compatible - within the accuracy of our data 
- with that found in the case of the distribution of the smallest eigenvalue, 
namely chUE in the large- volume region and chSE in the small- volume region. 

2 We thank T. Wettig for pointing this out to us. 



«ai« 

™qcd m n 

which, in our case with Nf = 1 reduces to 



(32) 




(33) 
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Figure 6: Distribution of the smallest eigenvalue -P(A m ; n ) in the v — 
topological sector for different lattices and values of j3, for the FP and Neu- 
berger's Dirac operator. The full and dashed curves represent the theoretical 
predictions of chUE and chSE, respectively. 

Number variance. A third test was accomplished by checking the number 
variance £ 2 (So, S). We report the results (for the unfolded variable z) in Fig. 
IH1 for the same collection of cases considered in Fig. |H1 and In the large 
volume case we see agreement with the theoretical prediction of chUE 
up to a value of z ~ 2. This outcome is also consistent with the results for 
p s (z), where the lattice data agree with chUE up to the same value of z, 
corresponding to the support of the distribution of the smallest eigenvalue 
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Figure 7: Microscopic spectral density p s (z) in the v — topological sector 
for different lattices and values of /3, for the FP and Neuberger's Dirac op- 
erator. The full and dashed curves represent the theoretical predictions of 
chUE and chSE, respectively. 



(see Fig. Inland [3 upper part). The existence of an upper bound of z, z max , 
for the applicability of chRMT is well-known in the literature (see for example 
48J) and an estimate for z max is available in d — 4: z max ~ f% {La) 2 ] 
such an estimate (related to the massless modes in d = 4) is absent in d — 2. 

In the small volume case we observe instead disagreement with the ana- 
lytical prediction of chSE ^jg^ already for very small values of z, even if the 
overall shape seems to resemble the symplectic behavior. 
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Figure 8: Number variance as a function of the unfolded variable z: (a) 
large volume region, L = 16, (5 = 2 and Neuberger's action (circles), L = 24, 
(3 = 4 and FP action (boxes); (b) small volume region, L = 16, (3 = 4 for 
Neuberger's (triangles) and the FP action (diamonds). The full and dashed 
curves represent the theoretical predictions of chUE and chSE, respectively. 

4.2 Topological sector: Nf = and v = 1 

As discussed earlier ( Section l2~3j) . in the case of GW actions a definition of the 
topological charge corresponding to the parameter v of chRMT is available: 
this is the index of the lattice Dirac operator. 

In Fig. Eland El we present our results for the distribution of the smallest 
eigenvalue and microscopic spectral density respectively, for configurations 
with index(i7) = l and compare with the corresponding predictions of chUE 
[oTH 121] and chSE [!J7J H3 EH] • For these configurations we use the identical 
unfolding as in the previous section for v = configurations (for correspond- 
ing L and 0) and, in particular, the same value of S. This procedure is 
consistent with the expectations from chRMT that S should not depend 
on the topological sector, since all such dependence is incorporated in the 
functional behavior of the distributions. We observe for P(A min ) and p s {z) 
a universal behavior, again consistent with the observations of the previous 
section. 
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Figure 9: Distribution P(z min ) of the smallest eigenvalue in the unfolded 
variable z in the v — 1 topological sector for different lattices and values of 
P, for the FP and Neuberger's Dirac operator. The dotted and dashed curves 
represent the theoretical predictions of the chRMT ensemble for the values 
v — and v — 1, respectively, with chUE for (a, b) and chSE for (c, d), cf. 
the discussion in the text. 
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Figure 10: Microscopic spectral density p s (z) in the v — 1 topological sector 
for different lattices and values of (3, for the FP and Neuberger's Dirac oper- 
ator. The dashed curves represent the theoretical predictions of the chRMT 
ensemble for v — 1, with chUE for (a, b) and chSE for (c, d). 



4.3 Dynamical fermions 

Consistent with our unquenching procedure, we define the distributions for 
the setup with dynamical fermions by weighting the entries of the quenched 
distributions with the fermion determinant. For example, in the case of the 
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microscopic spectral density we take: 



Ps(z) 



EcPs(C,z) detV(C) 
J2 c detV(C) 



(34) 



where p s (C,z) is the spectral distribution for an individual configuration 
C; the quenched setup is recovered replacing detV (C) with unity. This 
"unquenching" is formally correct; it is not optimal from the point of view of 
statistical reliability and one may expect large fluctuations of intermediate 
results. However, in our d = 2 studies (cf. [131 CHI) we found stable results 
for the given statistics. 

In Fig. flTI and fT2l we show the results for P(A min ) and p s (z) respectively 
in the case of configurations with index({7)=0. We compare these with the 
predictions of chUE [EIIl EH and chSE [SB 12 EE] in the case N f = 1, v = 0. 
In Fig. El we see that higher statistics for dynamical fermions would be 
desirable in particular for p s (z). 

The value of £ used for the definition of the unfolded variable z was found 
by a best-fit procedure in the distribution of the smallest eigenvalues like in 
Section l4~Tl In the large- volume case we find, consistent with the results of 
the quenched setup, agreement with chUE up to z ~ 2. 

The values of £ found here provide an estimate of the fermion condensate 
in the infinite volume limit. In Table El we summarize these results for this 
quantity. These should also be compared with the direct determination in 
the finite volume of Section EH"! in Table El Typically, the statistical error are 
smaller for the values obtained by the RMT techniques. 

In Fig. El we display these estimates for the infinite volume condensate. 
The results obtained in the large volume region (using the chUE distribu- 
tions) seem to approach the continuum with a linear dependence in the lattice 
spacing; those obtained in the small volume region appear more scattered. 

Another check [HI] of the chRMT predictions related to the microscopic 
spectral distribution involves the quantity [SU] 




(35) 



Inserting the predictions from chRMT for p s (z), one obtains [2*5] : 



8(** + f + JV/-1) 



(36) 
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Figure 11: Distribution P(z min ) of the smallest eigenvalue in the unfolded 
variable z for the unquenched (Nf = 1) data in the v — topological sector, 
for the FP and Neuberger's Dirac operator. The dotted and dashed curves 
represent the theoretical prediction of the best fitting chRMT ensemble (cf. 
Fig. El) in the trivial sector for Nf = and Nf = 1, respectively. 

with b = 2 and 4 for chUE and chSE, respectively. 

In Tables H] and El we compare our values for S2 with (|36J) in the large 
and the small (physical) volume region. We find better agreement for large 
volumes and for higher topological sectors. 
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Figure 12: Microscopic spectral density p s (z) for the unquenched (Nf = 
1) data in the trivial topological sector for the FP and Neuberger's Dirac 
operator. The dashed curves represent the theoretical predictions of the best 
fitting chRMT ensemble (cf. Fig. EJ) for Nf = and Nf = 1, respectively. 



5 Discussion and conclusions 

Both, the FP action and Neuberger's Dirac operator, guarantee the restora- 
tion of the main features of chiral symmetry in the continuum limit. We 
checked this through the analysis of their spectrum, in particular the zero 
modes and their chirality. The GW relation allows one to define an index for 
the Dirac operator which is useful to define a topological charge on the lat- 
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Figure 13: The ratio of E yf[3 (in units of the continuum value — (V>V)cont/e = 
c, cf. flUH) )) obtained from the distribution of the smallest eigenvalue for the 
quenched (left) and unquenched (right) setup as a function of ea = l/y/fi. 
Full (open) symbols denote results of the fit in the large (small) physical 
volume region, with unitary (symplectic) prediction; circles denote the FP 
action, boxes Neuberger's operator. 
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Table 3: Values of the fermion condensate in the infinite volume obtained 
from comparison of lattice data with chRMT. 



tice. The additional property of the FP action to be perfect at the classical 
level (which is not the case for Neuberger's operator) makes this definition 
particularly reliable, since it is related to a renormalization group guided 
procedure of interpolation of the gauge configuration. This results, as we 
observe, in an excellent agreement - good even at low (3 - with the geometric 
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Table 4: Our values for S2 in the large-volume region in the quenched and 
unquenched setup: L = 24 and (3 = 4 (FP action), L = 16 and (3 = 2 
(Neuberger's Dirac operator); chUE denotes the prediction (|3l)j) for 6 = 2. 
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Table 5: The values 5*2 in the small-volume region in the quenched and 
unquenched setup: L — 16 and (3 = 4; chSE denotes the prediction (|3^|) for 
b = 4. 



definition, which in this context is the most close-fitting to the continuum. 
The agreement is poor in the case of Neuberger's operator. 

At the practical level, Neuberger's approach seems to be most suitable 
when strictness of the chiral properties is required, as in the case of the 
calculation of the fermion condensate by subtraction. Parametrization un- 
certainties of the FP action seem to be a problem in this case, for small or 
even vanishing fermion mass. On the other hand, some experience ^Ej in- 
dicates that e.g. for dispersion relations, even the approximate perfectness 
of the FP action is a great advantage, while Neuberger's operator introduces 
discretization effects as large as for the Wilson action. Of course, the present 
parametrization of the FP action could be improved (e.g. by sampling ac- 
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cording to the compact gauge action in the determination of the FP fermion 
action); there may be even the chance to reduce the effective number of pa- 
rameters jSI]. However, in view of the four- dimensional environment, where 
only few operators can be included in the action for practical reasons, the 
parametrization effects are likely to be a serious problem, at least for small 
fermion mass. 

We also studied the statistical fluctuations of the spectra of the Dirac 
operators V Fp and V Ne • In all cases considered here (different values of 
(3 and lattice size) V Fp and £>Ne display very similar statistical properties. 
Concerning comparison with chRMT, the results show an unexpected com- 
plementarity between the symplectic and the unitary universal behavior, the 
latter coming into play for larger physical volumes. However, several aspects 
lead us to conclude that the symplectic behavior is an artifact of the small 
volume where RMT breaks down. The spectrum is not doubly degenerate 
as would be expected in the symplectic framework; it just looks like one half 
of a chSE spectrum. The number variance shows deviation from chSE al- 
ready from very small values of z. The distribution shape depends on the 
physical size of the volume. This is in agreement with the expected range 
of applicability of RMT, l/(am Schwinger ) <C L as discussed in Section H~T1 
We conclude that RMT necessitates randomness at a lattice size L much 
larger than the correlation length of the system. In the (in physical units) 
large volume region chUE describes our data up to z max ~ 2 (for the trivial 
sector). This distribution then should govern the continuum limit. 

RMT provides a powerful tool to extract the infinite volume chiral con- 
densate from finite (and hopefully small) lattices in the case of dynamical 
fermions. Our findings indicate, that some caution is in order. We argue that 
the physical volume should be not too small, since chRMT interpretation may 
produce misleading results in this situation. 

Concerning the condensate as determined from the chUE fits in the large 
volume domain, we observe linear cut-off effects for both actions (Fig. [T3*j) . 
If this is not an artifact of this approach, the natural explanation would be, 
that one still has to improve the field operators. This is also true for the FP 
action, since perfectness applies to spectral quantities and its extension to 
other observables passes through the implementation of the RG improvement 
to the operator (in this case ipip) as well. Also, we did not attempt to 
introduce a renormalized coupling constant. 

We find that studying the spectrum of GW operators is extremely helpful 
to identify chirality properties. With chRMT one learns how to disentangle 
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universal properties from the dynamical parameters like the fermion conden- 
sate. Our FP action and Neuberger's operator lead to similar results in this 
context. 
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